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\ 1. Introduction 

With three operational fluorescence sites out of four and more than 1000 
\^ ' active Cherenkov detectors on the ground at the time of writing, the Pierre 

Auger Observatory is nearing completion and has started accumulating 
data at a regularly increasing pace. In spite of the still small statistics 
Q_i' available, a lot of progress has been made in the understanding and fine- 

q . tuning of the detector, which has resulted in the development of reliable 

S-j \ analysis methods and of the release of its first scientific results concern- 

ing the main issues in ultra-high energy (UHE) cosmic ray physics. The 
spectrum of UHE cosmic rays observed by Auger is presented in 1 as an 
illustration of the power of Auger's hybrid detection, and the present con- 
■ tribution will focus on the results obtained in the context of anisotropy 

^ \ searches and composition studies. 

2. The arrival direction of UHECR: anisotropy studies with 
the Auger Observatory 

Anisotropics in the flux of UHE cosmic rays may appear in different energy 
ranges and angular scales, depending on the nature, distance and extension 
of the source. Cosmic rays around an EeV are thought to be of galactic 
origin, and the region of the Galactic Center and the Galactic Plane are 
key targets for anisotropy searches performed with Auger data. At higher 
energies one rather expects UHE cosmic rays to come from extra-galactic 
sources; a search for directional excesses of cosmic rays could then reveal 
a correlation with some (un)known astrophysical objects or even exotic 
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sources. 

The anisotropy studies performed by Auger are based on all surface 
detector (SD) events (plus some hybrids) with zenith angle 9 < 60° that 
pass the quality cut T5, which requires that the station with the highest 
signal be surrounded by a hexagon of working stations, ensuring a good 
reconstruction of the event. The energy of the events is determined using 
the constant intensity cut method and calibrating the parameter to the 
energy obtained from the florescence detector (FD) as described in. 1 

2.1. Angular resolution and coverage maps 

To detect an excess of events coming from a particular region of the sky, 
one has to compare the number of events observed in that region with the 
corresponding coverage map, that is, the background number of events ex- 
pected from an isotropic flux of cosmic rays in the same exposure conditions. 
This procedure requires accurate knowledge of the detector properties, and 
in particular of its angular accuracy and of the exposure dependance in 
time, energy and solid angle for each point of the sky. For a more detailed 
discussion of these issues, we refer the reader to. 2 ' 3 

The angular resolution AR for the SD is defined as the angular radius 
that would contain 68% of the showers coming from a point source; it is 
determined from the zenith (9) and azimuth (<f) uncertainties obtained from 
the geometrical reconstruction on an event by event basis, 

AR = 1.5 y/[a 2 tfi) + sin 2 {9) ct 2 (0)] /2 (1) 

where a 2 is for the variance. The AR is driven by the accuracy on the 
measurement of the arrival time of the shower front in each station. The 
variance on the arrival time T\ of the first particle is parameterized accord- 
ing to the time variance model described in, 2 ' 3 which was validated using 
data from the so-called "doublets" (pairs of tanks separated by 11 m). 

To build the coverage maps, one has to consider all possible modula- 
tions and inhomogeneities in the exposure of different regions of the sky. 
Besides the obvious effects due to the rotation of the Earth and the limited 
field of view of the detector, other modulations are induced by the contin- 
uously growing size of the array, by temporary failures in some detectors, 
and by temperature and pressure variations (which affect both the shower 
development in the atmosphere and the response of the electronics) . Two 
different techniques have been used to estimate the SD coverage maps: 4 

• the semi-analytic method consists in an analytical fit to the 9 
distribution of the events in the relevant energy range, convoluted 
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with an acceptance factor which accounts for the time evolution of 
the detector (according to the trigger activity), assuming a uniform 
response in azimuth (which is valid for showers up to 60°). 
• the shuffling technique takes the average of many fake data sets 
generated by shuffling the observed events in such a way that the 
arrival times are exchanged and the azimuths are drawn uniformly. 
This shuffling also preserves the 9 distribution of the events. It 
might partially absorb an intrinsic large-scale anisotropy present 
in the cosmic ray flux, but this drawback can be avoided using 
independent shufflings in (day x hour). 

The expected number of events in a given pixel of the sky is obtained 
by integrating the coverage map in a given window, while the signal is 
determined by applying the same filtering to the event map. A significance 
map is then generated by comparing the signal in each pixel respect to the 
expected background, according to the Li & Ma procedure. 5 

2.2. Anisotropy studies around the galactic center 

The region of the Galactic Center (GC, located at the equatorial coordinates 
(a, 5) = (266.3°, -29.0°)) and the Galactic Plane (GP) are particularly 
attractive targets for cosmic ray anisotropy studies around EeV energies. 
Two cosmic ray experiments, AGASA and SUGAR, have already claimed 
significant excesses in the flux of UHECR in that region. AGASA 6 reported 
a 4.5cr excess of CR with energies in the range 10 18 — 10 184 eV in a 20° 
radius region centered at (a, S) ~ (280°, —17°) (it is worth noting however 
that the GC itself lies outside of the AGASA field of view). Subsequent 
searches near this region using old SUGAR data 7 failed to confirm that 
result but found a 2.9a excess flux of CR in the energy range 10 179 — 10 18 5 
eV in a 5.5° window centered at (a, 5) ~ (274°, —22°). Recent observations 
by HESS of a TeV 7 ray source in that region 8 and of diffuse 7-ray emission 
from the central 200 pc of the GP 9 have provided additional hints towards 
the presence of powerful CR accelerators in the Galaxy. In that context, 
several models that predict a detectable flux of neutrons in the EeV range 
(whose decay length is about the distance from the GC to the Earth) have 
also been proposed. 

With the GC well in the field of view and an angular resolution which is 
much better than previous CR experiments, the Pierre Auger Observatory 
is well suited to look for UHECR anisotropies coming from that region. A 
total of 79265 SD events and 3934 hybrid events have been used, which 
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Fig. 1. Left: significance map of CR overdensities in the region of the Galactic Center 
in the energy range 10 17 ' 9 — 10 18 ' 5 eV, showing the Galactic Center (cross), the Galactic 
Plane (solid line), the regions of excesses of AGASA and SUGAR (circles), and the 
AG AS A field of view limit (dashed line). The event map was smoothed with a top-hat 
5.5° window. Bight: corresponding histogram of overdensities computed on a grid of 3° 
spacing, compared to the average isotropic expectations points (with 2<r bars), (from 10 ) 



corresponds to the data collected between January 2004 and March 2006 
satisfying the T5 quality cut 1 and with 9 < 60°, 10 179 eV < E < 10 185 eV; 
it represents respectively more than four and ten times the sample that 
AGASA and SUGAR used in this context. 

Significance maps were built using different filterings of the data to ac- 
count for the angular size of the excesses reported by AGASA and SUGAR 
in their respective energy range. An example of such a map is shown in Fig. 
1 together with the corresponding overdensity distribution. Several tests 
were also performed with modified energy windows to account for a possi- 
ble energy shift due to differences in the calibration of the experiments. In 
all cases, Auger data have been found compatible with isotropy, therefore 
not confirming the results from previous experiments. Even in the worst 
case of a source emitting nucleons and embedded in a background made 
of heavier nuclei, to which Auger is more sensitive in the relevant energy 
range, a significant excess ( 5.2a) would be expected, in contradiction with 
current observations. Details of the analysis can be found in. 10 

Data from the Auger Observatory were also used to search for a point 
source in the direction of the GC itself at the scale of Auger's own an- 
gular resolution. In the energy range I0 179 - I0 185 eV, and applying a 
1.5° Gaussian filter to account for the pointing accuracy of the SD, we 
obtain 53.8 observed events against 45.8 expected. This allows to put a 
95% C.L. upper bound on the number of events coming from the source 



February 5, 2008 15:10 



WSPC - Proceedings Trim Size: 9in x 6in VVE'AugerResults 



5 

of n® 5 = 18.5. Assuming that, in the energy range considered, both 
the source and the bulk CR spectrum have similar spectral indexes and 
that the emitted CR are proton-like, and taking a differential spectrum 
$Cr{E) — £ 30 (E I EeV)~ 3 EeV~ 1 km~ 2 yr~ 1 sr~ 1 , where £ parameterizes 
the uncertainties on the flux normalization, a 95% C.L. upper bound of 

$f < £ 0.08 km-Zyr- 1 (2) 

can be set on the source flux. This bound could however be about 30% 
higher if the CR composition at EeV were heavy, i.e. close to Iron. 

Finally, a scan for correlations of CR arrival directions with the Galac- 
tic Plane and Super-Galactic Plane have also been made in two different 
windows of energy (1 EeV < E < 5 EeV and E > 5 EeV), yielding again 
negative results (although with a smaller dataset). 11 

2.3. Other searches for localized excesses in the Auger sky 
maps 

The Auger data have also been used to perform both blind searches and 
prescripted searches for localized excesses in other parts of the sky. 

In the case of blind searches, the distribution of significances is compared 
to those obtained from a large number of Monte Carlo isotropic simulations. 
Such searches were performed both for a 5° and a 15° angular scale and 
in two separate energy ranges, lEeV < E < 5EeV and E > 5EeV; all of 
them turned out to be compatible with isotropy. 12 

The Auger Collaboration had also released a list of prescribed targets 
with definite angular and energy windows, 13 with the associated signifi- 
cance probability level to attain in order to claim a positive signal. The 
prescription targets range from the Galactic Center to some nearby violent 
extragalactic objects; none of them has turned out to lead to a positive 
detection. As more data is streaming in, the catalogue of candidate targets 
that will be studied is expected to increase in the future. 

3. The nature of UHECR: composition studies with the 
Auger Observatory 

Thanks to its hybrid capabilities, the Auger Observatory can extract com- 
plementary information on the shower development parameters, that are 
ultimately related to the nature of the primary cosmic rays. If the discrimi- 
nation between different types of nuclei is complicated by the uncertainties 
in the hadronic models governing the interactions of the particles in the 
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Fig. 2. X max distribution of simulated 
and real candidate photon events. The 
blue point is the X ma x value and 
uncertainty for one event from the data. 
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Fig. 3. 95% C.L upper limit on the 
photon fraction in UHE cosmic rays 
obtained by Auger, compared to the 
results of Haverah Park, 15 HP, and 
AGASA, 16 Al and A2. 



shower at such high energies, several methods have already been proposed 
for the identification of photons and neutrinos, and are currently applied to 
the data of the Auger Observatory. The presence and the amount of pho- 
tons and neutrinos at such high energies would constitute a crucial probe 
for many exotic models of UHE cosmic ray production and could help locate 
candidate sources as they travel undeflected by the intergalactic magnetic 
fields. 

3.1. Upper limit on the UHE photon flux 

Unlike protons and nuclei, the development of photon showers are driven by 
electromagnetic (EM) interactions and do not suffer much from the uncer- 
tainties in hadronic interactions. Photon showers are expected to contain 
fewer and less energetic secondary muons, as a result of the smallness of 
the photon radiation length respect to its mean free path for photo-nuclear 
interactions and direct muon production. Their development is also delayed 
due to the small multiplicity in EM interactions and to the LPM effect, 14 
which reduces the bremsstrahlung and pair production cross-sections at 
energies above 10 EeV. These considerations allowed several ground array 
experiments to set upper limits on the flux of UHE photons on basis of 
studies of the rate of vertical to inclined showers (in Haverah Park experi- 
ment 15 ) and of the muon content of the showers at ground (in AGASA 16 ). 

Taking profit of its hybrid design, Auger has set up a different method 
to identify photon primaries in the flux of UHECR. It is based on the 
direct observation of the longitudinal profile of the shower development 
in the atmosphere by the FD, and uses as a discriminating variable the 
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atmospheric depth of the shower maximum, X max (the estimated average 
difference in X max between photons and hadrons is about 200gr/cm 2 ). 

The data set used for this analysis corresponds to the hybrid events 
(i.e. those observed by one or more FD telescope and by at least one SD 
station, which ensures a better angular accuracy and smaller uncertainty 
in the reconstruction of X max ) with a reconstructed energy E > 10 19 eV, 
registered between January 2004 and February 2006. During that period 
two of the four Auger eyes were active (for a total of 12 FD telescopes) and 
the number of deployed SD stations grew from 150 to 950. 

A series of cuts were applied to the data that guarantee the quality of 
the hybrid geometry and of the fit to the shower longitudinal profile, which 
takes into account the local amtospheric conditions (see the detail in 17 ). One 
important condition is to have the X max of the shower inside the field of 
view of the telescopes. To minimize the bias that this condition introduces 
against photon primaries in the detector acceptance, additional energy- 
dependant cuts are applied both on the zenith angle and the maximum 
distance of telescope to shower impact point in order to eliminate nearly- 
vertical and distant events. 

For each of the 29 events that survived all the cuts, 100 photon show- 
ers were simulated in the same energy and arrival direction conditions and 
the resulting expected distribution of X max was compared to the observed 
X m ax of the event. An example is shown in Fig. 2, together with the distri- 
bution of the X max from the whole selected dataset. For all 29 events, the 
observed X. max is well below the average value expected for photons. Taking 
systematic uncertainties on the X max determination and the photon shower 
simulations into account, the available statistics allows to put an upper limit 
on the photon fraction of 16% at 95% C.L, which is shown in Fig. 3 together 
with previous results and some predictions from non-accelerator models. 

3.2. Inclined air showers and the detection of neutrinos 

The use of Cherenkov water tanks for the SD allows the Pierre Auger Ob- 
servatory to detect showers with zenith angles up to 90° (and even more). 18 
The range of inclined showers, 60° < 9 < 90°, contributes half the total 
solid angle of the detector and about 25% of its geometrical acceptance, 
thereby significantly increasing the field of view of the detector and the SD 
statistics. Such events are indeed seen by Auger both in the SD and the 
FD; some of them may be quite spectacular, with very extended footprints 
involving tens of tanks, as illustrated in Fig. 4. Dedicated selection pro- 
cedures and reconstruction methods are being developed in Auger to deal 
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Fig. 4. Example of a near horizontal air shower as seen by the SD; the shower triggered 
31 tanks and extends on about 30 km at ground. In the upper left corner, the best fitting 
simulated muon map corresponding to the reconstructed zenith and azimuth angles. 

with the distinctive features of those showers. 

The distance between the first interaction point (normally in the first 
few 100 g cm -2 ) and the detector position is much larger than in the vertical 
case, the atmospheric depth ranging from 1740 g cm' 2 at 60° till 31 000 
g cm~ 2 at 90°). As a result, the EM component of the shower dies out 
long before reaching the ground, and the only particles recorded in the SD 
are energetic muons (typically of 10-1000 GeV) accompanied by an EM 
halo which is constantly regenerated by muon decay, brehmsstrahlung and 
pair production. Those muons arrive at ground in a thin front with small 
curvature, resulting in short FADC pulses in the tanks, as shown in Fig. 5 
(right). Their trajectories are long enough to be affected by the geomagnetic 
field, which leads to a separation between positive and negative muons and 
a further elongation of the projected footprint on the ground. 

The reconstruction of inclined showers is based on the search for the 
best fit to the pattern of signals at ground performed with averaged maps of 
muon densities obtained from simulations. The relation between the muon 
density and the energy depends on the nature of the primary cosmic ray, 
and is established on basis of Monte Carlo simulations which suffer from 
hadronic interactions uncertainties at high energies. In this context, hybrid 
inclined events reconstructed by both the SD and the FD will play an im- 
portant role in primary composition studies, since they allow independant 
measurements of the EM and muonic components of the shower. 19 

Inclined showers also constitute the bulk of events from which a signal of 
UHE cosmic neutrino could be extracted. Due to their small cross-section, 
neutrinos can penetrate deeply in the atmosphere and initiate showers at all 
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Fig. 5. FADC traces of a young (left) and old (right) shower. The signal from a young 
shower gets smaller and more extended as the distance to the core increases, while old 
showers have short traces at all distances. 

possible depths, unlike nuclei or photons. In particular, showers originating 
less than w 2000 g cm" 2 away from the detector will reach it before their 
EM component attenuates completely. Selection criteria will thus require 
the presence of signals corresponding to a young shower, and in particular of 
stations with extended traces that reflect the large curvature of the shower 
front and the presence of an EM component (see Fig. 5). 

Up-going tau neutrinos that skim the Earth just below the horizon could 
also be detected as they are likely to interact in the ground and produce 
a tau which may emerge from Earth and initiate an observable air shower, 
provided it decays close enough to the SD. Preliminary studies provided a 
proof of principle for the detection of such neutrinos in the energy range 

10 17 10 19 eV 20 and, although a careful study of systematic uncertainties 

is necessary to infer with a reasonable precision the energy of the incident 
v T primary, this method seems the most promising in terms of acceptance, 
which is a crucial matter when dealing with event rates as small as ~ 1 per 
year. Studies are currently ongoing both in the down-going and upgoing 
ranges to define and optimize the selection criteria, and the search for UHE 
neutrinos in the Auger data has started. 

4. Conclusions 

The Southern Auger Observatory, expected to be complete in 2007, has de- 
livered its first science results on the UHE cosmic ray spectrum, anisotropy 
searches and composition studies. In particular, the region of the Galactic 
Center has been studied with a precision never attained before, yielding no 
hint of anisotropics. The absence of evidence for a point-source near the GC 
excludes several scenarios of neutron sources recently proposed. The upper 
limit on the photon fraction above 10 EeV, derived for the first time from 
a direct observation of the shower maximum, confirms and improves pre- 
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vious limits from ground arrays. Finally, the inclined shower data sample 
will soon contribute to enlarge the field of view of the detector and increase 
its statistics; it might also reveal the first cosmic neutrino ever observed at 
ultra-high energies. 
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